Cooling in the single-photon strong-coupling regime of cavity optomechanics 
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In this paper we discuss how red-sideband cooling is modified in the single-photon strong-coupling regime 
of cavity optomechanics where the radiation pressure of a single photon displaces the mechanical oscillator 
by more than its zero-point uncertainty. Using Fermi's Golden rule we calculate the transition rates induced 
by the optical drive without linearizing the optomechanical interaction. In the resolved-sideband limit we find 
multiple-phonon cooling resonances for strong single-photon coupling that lead to non-thermal steady states 
including the possibility of phonon anti-bunching. Our study generalizes the standard linear cooling theory. 
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Introduction. In optomechanical systems mechanical de- 
grees of freedom are coupled to modes of the electromagnetic 
field inside an optical or microwave resonator |[ll|2l- Possible 
applications include ultra-sensitive sensing of masses, forces 
and electromagnetic fields |3], transducing quantum informa- 
tion between different parts of quantum networks and ex- 
ploring decoherence at larger mass and length scales fs*]. 

For these applications it is very important to minimize the 
influence of thermal fluctuations. This is why a large part of 
current experimental efforts is directed at cooling the mechan- 
ical degrees of freedom. Recently, the quantum ground state 
of mechanical motion was achieved for mesoscopic oscilla- 
tors [6-8J and the zero-point motion was detected by observ- 
ing an asymmetry between phonon absorption and emission 
rates, originating from the fact that an oscillator in the quan- 
tum ground state cannot emit but only absorb energy [9, 10]. 

In most optomechanical setups the position of the mechan- 
ical oscillator linearly modulates the cavity frequency. While 
the optomechanical coupling on the single-photon level is usu- 
ally much smaller than the cavity linewidth, coupling between 
the fluctuations of the cavity field and the mechanical position 
can be made appreciable using a strong optical drive. For this 
situation a quantum theory of red-sideband cooling has been 
developed in Refs. lilTlfl^. 



Several experiments |17|, |9|, [IJ, [141] are currently approach- 
ing the regime where the presence of a single photon displaces 
the mechanical oscillator by more than its zero-point uncer- 
tainty. Going beyond early work ifTsi [T6ll novel effects in this 
regime have recently been predicted, including mechanically- 
induced cavity resonances [17, 18], multiple mechanicalside- 
bands [17], photon anti-bunching [18], non-Gaussian 1 1 7i [l9ll 
or non-classical [20[ mechanical steady-states, and scattering 
ll2Tll of and interferometry [22[ with single photons. However, 
cooling of the mechanical oscillator in the regime of nonlinear 
strong coupling has not been discussed in the literature. 

In this article we study how the weak-coupling cooling the- 
ory is modified in the single-photon strong-coupling regime. 
Using Fermi's Golden rule we calculate the transition rates 
caused by the coupling to the optical field without lineariz- 
ing the optomechanical interaction. In the resolved-sideband 
limit we find cooling resonances if the cavity is driven on one 
of the several mechanical sidebands. In contrast to the weak- 
coupling regime the phonon transition rates do not obey de- 



tailed balance. We find steady states with non-thermal phonon 
number statistics including phonon anti-bunching. 

Our study generahzes the standard theory of red-sideband 
cooling lllll[l2[l to the regime of strong optomechanical cou- 
pling. In the literature nonlinear cooling has been discussed 
in the context of trapped ions outside the Lamb-Dicke regime 
and of optomechancial systems where the cavity is cou- 
pled to the position squared of the oscillator ll24ll . 

Model. We consider the standard model of optomechanical 
systems where the position of a mechanical oscillator, x ~ 
2^zpf(& + b''), is parametrically coupled to an optical cavity 
mode a. The Hamiltonian reads 



Hq ~ hjjjia) a + hbJMb'^b + hga^ a{b + 6^) 



(1) 



where ujn is the resonator frequency, lom the mechanical fre- 
quency, and g = w^xzpf is the optomechanical coupling. 
2^zPF = \/ f^l {'^Mujm) is the zero-point uncertainty, M the 
mass of the mechanical oscillator, and uj'r, = the deriva- 
tive of the resonator energy with respect to the oscillator po- 
sition X. a and h are bosonic annihilation operators for the 
cavity mode and the mechanical oscillator, respectively. 

In order to include drive and decay we use standard input- 
output theory f23\ . In a frame rotating at the frequency of the 
optical drive, the nonlinear quantum Langevin equations read 

a = +iAa — —a ~ ig (^b^ + bj a + V^coin (2) 



b = —iojMb — —b — igcva + ^J^b^ 



(3) 



where A = oj^ — loji is the detuning between laser wl and 
resonator frequency ojn, and 7 and k are the mechanical and 
cavity damping rates. The cavity input Sin = Sin + ^ is the 
sum of a coherent amplitude a^^ and a vacuum noise operator 
i satisfying m^^{t')) - 5{t - t') and {CHmt')) = 0. 
Finally, we assume that the mechanical bath is Markovian and 
has a temperature T, i.e. (Sin(<)6j„(*')) = ("th + - t') 
and {bl^it)bin{t')) = nti,S{t~t') with n^^^ = griu^M/kBT ^ 
Multiple cooling resonances. We are interested in the influ- 
ence of a weakly driven, but strongly coupled optical field on 
the mechanical oscillator. We calculate the mechanical cool- 
ing and amplification rates induced by the optical drive using 
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FIG. 1: (Color online) Multiple cooling resonances. Steady-state 
phonon number n = {o'b) (a) for q/lum ~ 0.5 (red solid), g/wM ~ 
0.25 (green dashed) and g/uM = 0.1 (blue dash-dotted) as well as 
/k = 4 and (b) for ujm /hi = 0.5 (red solid), uim / hi = 2 (green 
dashed) and cjAf/ft = 4 (blue dash-dotted) as well as gjiou ~ 0.5. 
Results from the set of rate equations l|4) are shown as lines and those 
from the quantum master equation ([8} as dots. The other parameters 
are liJmH ~ 1000, nth = 1> and 0,1 k — 0.2. 



Fermi's Golden rule and write down a set of rate equations for 
the mechanical oscillator 

Pn = -^nth{n + 1)P„ - 7(nth + l)nP„ 
+ jnthnPn-i + j{nth + !)("- + l)-Pn+i 



where P„ is the probability of the oscillator to be in the state 
with n phonons. The terms in the first two lines are due to the 
coupling of the mechanical oscillator to its thermal bath with 
rate 7 and thermal phonon number nth. The sums in the last 
line are the terms caused by the coupling to the cavity field. 

In the frame rotating at the drive frequency ljl the drive is 
described by the Hamiltonian H[ = fi,f2(a + a^) with fl = 
\/K|ain|- We seek the transition rates r„_j.m/n from the state 
\n) with n phonons to the state |m) with m phonons induced 
by the optical drive. We work to second order in H, to obtain 
a Fermi Golden rule result [25,1 valid for O <C k where cavity 
states with more than one photon can be neglected. This gives 

r„^rn#« =-^^rfri^dT2(n|(z|i7((Ti)|m)(m|i7((T2)|z)|n) 

t Jo Jo J-00 J-00 

g-(K/2-iA)(Ti-Si)g-(K/2 + iA)(T2-S2) 



m) 



(5) 



where A = A + g'^/wM, |«) is the vacuum state of the optical 
bath, and we used the solution to Eq. (|2]i in the absence of an 
optical drive, as derived in the appendix. 

Using a resolution of unity we rewrite the matrix element 

(n|e-*(^i)e-*('^i)|m) = J2k{Me^^^'^\'^){k\e'^^"'^\m). For 
large mechanical quality factors we only need to consider the 

free mechanical evolution {n\e^^^^'>\k) = e'("-'=)"»^"'iZ„,fe 
where we have evaluated the matrix elements to be Z„ k = 



min(n,fc)! j^(\n — k\) 



x(n,/c)l 



with the associated Laguerre polynomials i„ 



n(n,fc) 

x) 



268 and 



the coupling strength A = g/ujM- Finally, we obtain the rates 
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^71 , k , k 



y 

^K/2-i[(n-A:)a;A/ + A] 



(6) 



In the resolved-sideband limit um ^ only terms with k - 
n — / > contribute significantly for detunings A w —Iujm 



^0272 72 
(«;/2)2 + (luJM + A)2 



(7) 



Eqs. (|6]l and O are our main result from which we obtain a 
clear physical picture how an incident photon is inelastically 
scattered off the cavity and changes the state of the mechani- 
cal oscillator from n to to phonons. In the resolved-sideband 
limit lum ^ K and for a drive detuned by A = —Ilom, the 
process of destroying I phonons when creating the cavity pho- 
ton is enhanced by the cavity susceptibility. The amplitude for 
this process is proportional to the matrix element Z„ — 
J dx (p'^_i{x — xo)(pn{x) where (p„i{x) are the eigenfunctions 
of the simple harmonic oscillator and xo = —2xzPFgl^M is 
the displacement caused by a single photon. That means that 
the matrix element is given by an overlap integral between dis- 
placed harmonic oscillator wave functions in accordance with 
the Franck-Condon principle. As the photon leaves the cavity 
it induces a transition in the mechanical oscillator from n — I 
to TO phonons. This process is not resonantly enhanced as 
the photon decays into the continuum of modes in free space. 
This is why its amplitude is just given by the matrix element 
Zn-i,m which is a function of the ratio = g"^ i.e. the 
strength of the optomechanical interaction g relative to the fre- 
quency of the mechanical oscillator lom, but does not depend 
on the drive detuning A. The photon in the output field has 
an energy ftw^ + (n — m)hujM, i-e. it carries away the en- 
ergy of 71 — TO phonons. In the non-resolved sideband limit 
^ K processes with different intermediate phonon num- 
ber n — k I contribute, and their amplitudes interfere ac- 
cording to Eq. 

In Fig.[T]we plot the steady-state phonon number n = {b^b) 
as a function of detuning A = ul — uj^ for different coupling 
strengths g/uiM (a) and sideband parameters uji\j / k (b). In (a) 
we observe that for weak drive and in the resolved-sideband 
limit ujm ^ K several cooling resonances appear when the 
detuning matches A — —Ilom with / integer For smaller 
sideband parameters (b) we notice that resonances merge. 

Validity of the rate equation approach. To investigate the 
validity of the set of rate equations (HI we solve numerically 
the quantum master equation 



Q 



H',g +KD[a\Q + -i{\ + nth)V[b]Q- 



-inth'D[b'^]Q 
(8) 

withH' = Hq+H'i where 2?[o]gi = dgd^ -{d^dg+ go'' d)/2 is 
the standard dissipator in Lindblad form and Hq is the Hamil- 
tonian ([T]i in the frame rotating at the drive frequency lul. 
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FIG. 2: (Color online) Validity of the rate equation approach and 
phonon anti-bunching. (a) Mean photon number {a' a) (blue crosses) 
and mean phonon number n = (b^b) (red solid) as well as (b) phonon 
number fluctuations F — {¥¥bb) / (b^b)'^ as a function of the drive 
strength Q for a detuning A/ujm — —2. Results from the set of 
rate equations © are shown as lines and those from the quantum 
master equation ([8} as dots and crosses. The other parameters are 
(jJAi/f^ = 4, g/oJM = 0.5, n/K = 0.2, anda;A//7 = 1000. 



In Fig. [U we plot the master equation ([8| results alongside 
with those from the set of rate equations (|4|i. The agreement 
between the two is excellent. The small deviations at detun- 
ings A w between the rate equation and master equation 
results stem from the rt-photon resonances at A = —ng^/ujM 
ITtIi . In their vicinity off-diagonal elements of the density ma- 
trix cannot be neglected and the rate equation approach fails. 

In Fig.|2](a) we plot the mean photon and phonon number, 
(d^a) and n — {¥h), as a function of the drive strength VI. 
Driving on the second red sideband, i.e. A /ujm = —2, we find 
that for small fl/n the phonon number n decreases quadrati- 
cally in the drive strength as expected. For large fl/n the 
set of rate equations ^ predict a finite minimal phonon num- 
ber determined by the optically induced rates However, in 
Fig.|2](a) we find that as the photon number {a^a) increases 
with stronger drive and the master equation result ^ starts 
to deviate from solution of the rate equations 

Non-thermal steady states. The rate equations (HI in general 
do not obey detailed balance and so their steady state will not 
be a thermal state. In Fig.[3]we plot the mean phonon number 
n = (b^b) and the number fluctuations F — (b^ bb) / b)'^ 
as a function of detuning A. We note that the cooling power 
n/nth of the various cooling resonances depends on the ther- 
mal phonon number nth. This is a consequence of the nonlin- 
ear dependence of the rates (|6]l on phonon number. Close to 
the cooling resonances the mechanical state clearly deviates 
from a thermal state whose number fluctuations are given by 
F ^ 2. We find both reduced F < 2 and enhanced number 
fluctuations F > 2. In Fig. [3] we also plot the mean phonon 
number n and number fluctuations as a function of the ther- 
mal phonon number nth for A = — wm and A = —2ujm- The 
mean phonon number n is a nonlinear function of the thermal 
phonon number nth and the fluctuations F can change from 
_F<2toF>2asa function of nth- In Fig-|2](b) we find 
that even phonon anti-bunching, i.e. F < 1, can occur. 

To understand this behavior let us look at the properties of 
the rate equations in more detail. For example, the resonant 
one-phonon cooling r„^„_i and amplification r„^„+i rates 
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FIG. 3: (Color online) Non-thermal steady-states. (Left) Phonon 
number n = (b^b) (top) and phonon number fluctuations F = 
(b^ b^ bb) / (b'b)'^ (bottom) as a function of detuning A for nth = 1 
(solid red) and nth = 10 (dashed blue). (Right) Phonon number 
n (top) and phonon number fluctuations F (bottom) as a function 
of the thermal phonon number nth for A — —lum (solid red) and 
A = —2ijJm (dashed blue). The other parameters are ljm ~ 4, 
g/ujM ~ 0.5, Q/k = 0.2, and ujm/"/ = 1000. Thin black lines 
show results of the linear model ( jilt and ( I12t for Q,X/k = 0.2. 



in the resolved-sideband limit ojm ^ k O read 



and 



r 



(n + l)ri_,o 



[£.W(A^)4°ii(A^)]^ 
(n + l)2 



(9) 



(10) 



In the special case of weak coupling A ^ 1 where one- 
phonon processes are most important, since higher-order pro- 
cesses are suppressed by a larger power of A, we obtain 



and 



r„- 



(^/2)2 + (A + wm)2 



(«/2)2 + (A - ujMr 



(11) 



(12) 



That means we recover the standard cooling theory iflll [T2I1 
where the rates are linear in n and n + 1, respectively, i.e. one 
can write them as r„_j.„_i = nFj^ and r„^„+i — {n + l)T^. 
In this case and simply renormalize the thermal mean 
phonon number n = (7nth + Ff) /{-y + Fj. — Ff ). Thus, the 
steady state is a thermal state with F = 2 for all detunings A 
and the mean phonon number n is linear in the thermal phonon 
number 7ith. In Fig.|3]we plot the mean phonon number n and 
the phonon number fluctuations F for the weak-coupling limit 
A <C 1. We find a single cooling resonance at A = —lum and 
F = 2, indicative of a thermal state. 
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In general, the normalized rates (|9]l and ( fTOl l differ from 
unity and multiple-phonon processes r„^m with |n — to| > 1 
are important. If only a few phonons remain, further insight 
can be obtained. For example, in Fig.|2](b) for 0/k w 2, only 
the states with zero, one and two phonons are important. The 
two-phonon cooling process r2^o reduces the occupation in 
the two-phonon state leading to phonon anti-bunching F < 1 
similar to the case studied in Ref. |24]. In contrast, in Fig. [3] 
for nth = 1 and A ~ — wm, we have r2_>i/(2ri^o) < 1> 
i.e. relative to a thermal state the one-phonon state is depleted 
faster than the two-phonon state which results in > 2. 

Detection. The multiple-phonon cooling and amplification 
processes lead to multiple mechanical sidebands in the opti- 
cal output spectrum lll7ll . This is not a proof of non-thermal 
states by itself, as it can also occur e.g. in the case of large me- 
chanical amplitude motion at weak optomechanical coupling 
ll27ll . The mean phonon number n and the phonon number 
fluctuations F can be obtained from the Wigner function of 
the mechanical oscillator. Schemes to reconstruct the Wigner 
function experimentally rely on back-action evasion li28ll . cou- 
pling to a two-level system || 29[ J s oll . or the time-dependence 
induced by short optical pulses Bill . 
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Appendix. To calculate the optically induced transition rates 
between different phonon Fock states (|5]l, we need an expres- 
sion for the optical field in the absence of an optical drive. The 
formal solution to Eq. ^ is 



aify — a/k / dTK{t,T)ain{T) 



(13) 



where A'(t,T) = e^^'^/^^^'^^^'^^^T [e^''^-''- 
and T is the time ordering operator. In the case of no optical 
drive, ain = 0, the operator identity a^a = holds, such that 
b{t') = e-*"^''(*'-*)S(t) for times \t' - i| < (7nth)"^ This 
means that for times (t — r) <C (7".th)^^, we can express the 
time ordered exponential above as e*^ {i-'^)/'^M gX(t)^-x(T) 
where X{t) — g[b{t) — W{t)]/LjjM- This follows from using 
the standard commutation relations for the bosonic operator b. 
Thus, in the limit 7nth ^ k, we find 



where A = + g'^ /ujm- This expression can also be derived 
using the polaron transform ITtIi . 
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